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Abstract 

The signatures of UHE proton propagation through CMB radiation are pair- 
production dip and GZK cutoff. The visible manifestations of these two spec- 
tral features are ankle, which is intrinsic part of the dip, beginning of GZK 
cutoff in the differential spectrum and -E1/2 in integral spectrum. Observed 
practically in all experiments since 1963, the ankle is usually interpreted as a 
feature caused by transition from galactic to extragalactic cosmic rays. Us- 
ing the mass composition measured by HiRes, Telescope Array and Auger 
detectors at energy (1 - 3) EeV, calculated anisotropy of galactic cosmic rays 
at these energies, and the elongation curves we strongly argue against the 
interpretation of the ankle given above. The transition must occur at lower 
energy, most probably at the second knee as the dip model predicts. The 
other prediction of the dip model, the shape of the dip, is well confirmed by 
HiRes, Telescope Array (TA), AGASA and Yakutsk detectors, and, after re- 
calibration of energies, by Auger detector. Predicted beginning of GZK cutoff 
and El 12 agree well with HiRes and TA data. However, directly measured 
mass composition remains a puzzle. While HiRes and TA detectors observe 
the proton-dominated mass composition, as required by the dip model, the 
data of Auger detector strongly evidence for nuclei mass composition be- 
coming progressively heavier at energy higher than 4 EeV and reaching Iron 
at energy about 35 EeV. The Auger-based scenario is consistent with an- 
other interpretation of the ankle at energy Ea ^ ^ EeV as transition from 
extragalactic protons to extragalactic nuclei. The heavy - nuclei dominance 
at higher energies may be provided by low-energy of acceleration for protons 
^max ^ ^ YjqY and rigidity-dependent E^^^ = ZE"^^^ for nuclei. The highest 
energy suppression may be explained as nuclei-destroying cutoff. 
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1. Introduction: spectrum features and signatures 

Do observed cosmic rays (CR) have galactic or extragalactic origin? 
Since 1950s this question divided theorists in two groups. The Russian group, 
presented most noticeably by V.L. Ginzburg [H and Syrovatsky defended 
the Galactic origin, while F. Hoyle jsj, G.R. Burbidge j3] and their collab- 
orators suggested the extragalactic origin. In the book [sj Ginzburg and 
Syrovatsky presented the strong arguments based on energetics, but hot dis- 
cussions between these two groups continued at all International Cosmic 
Ray Conferences (ICRC) until 1990s. V.L. Ginzburg suggested the rigor- 
ous test for extragalactic origin, predicting the low limit for gamma-ray flux 
from Small Magellanic Cloud (SMC). The gas density in this source is fairly 
well known, while the CR flux in extragalactic model is that measured in 
Milky Way (MW). Therefore the produced gamma-ray flux at energies 70 - 
100 MeV is well predicted. The calculated flux from pp — t- 7r° production 
gives the lower limit, because part of the detected flux can be produced by 
electrons and by the point-like sources. In 1993 EGRET (gI detected flux 
below the Ginzburg limit, which means that extragalactic CR flux is lower 
than that observed in MW. In 2008 Fermi/LAT 0] conflrmed this result with 
a high accuracy. The results of EGRET and FERMI/LAT are valid up to 1 
- 10 PeV. At higher (or much higher) energies the transition to extragalactic 
Cosmic Rays might be expected. 

The observed energy spectrum of Cosmic Rays (CR) has an approxi- 
mately power- law behavior for 11 orders of magnitude in energy with several 
features that can be linked with particles propagation and acceleration. 

The most prominent spectral feature is the knee in all-particle spectrum 
at energy 3-4 PeV, discovered flrst at the MSU (Moscow State University) 
array in 1958 j^. At the knee the spectrum E''^ steepens from 7 ^ 2.7 to 
7 ^ 3.1. This knee is provided by the light elements, protons and Helium, 
and is explained in the framework of the Standard Model (SM) for Galactic 
Cosmic Rays (GCR) by the maximum energy -Emax of acceleration in the 
Galactic Sources (Supernovae Remnants, SNR). In the case of the rigidity- 
dependent acceleration -Emax oc Z, where Z is charge number of a nuclei, the 
maximum acceleration energy is reached by Iron nuclei, and thus the Iron 
knee is predicted to be located at energy by factor 26 times higher than for 
proton knee, i.e. at energy -E^ax ~ (§0 ~ 100) PeV. Recently, the Iron knee 
was found indeed at energy 80 PeV in KASCADE-Grande experiment ^ in 
a good agreement with rigidity acceleration prediction. 
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The Standard Model (SM) for Galactic Cosmic Rays was first put for- 
ward by Ginzburg and Syrovatskii (sf and in the main features it was ac- 



complished in 1977 - 1978 by discovery of diffusive shock acceleration [10 
including amplification of magnetic field due to cosmic-rays streaming insta- 
bility (A.R. Bell 1978 in 10|). The SM is based on (i) supernova remnants 
as sources, (ii) SNR shock acceleration, and (iii) diffusive propagation of CR 
in the Galactic magnetic fields. 

Above the knee, at energy E'skn ~ (0.4 — 0.7) EeV, there is a faint feature 
in the spectrum |ll| called the second knee. It is seen in many experiments 
(for a review see [12|). This feature is often interpreted as a place of tran- 
sition from galactic to extragalactic CRs. However, for the last forty years 




PRIMARY ENERGY (ELECTRON VOLTS) 



Figure 1: The ankle as presented by J. Linsley in review |13| and taken from paper by 
T.K. Gaisser Q. 



the standard place for transition from galactic to extragalactic CRs is con- 
sidered at ankle, a very prominent spectral feature, observed first in 1960s by 
Vulcano Ranch detector at energy E°^^ ~ 10 EeV, and it was immediately 
interpreted by J. Linsley 15|] as transition between these two components of 
CRs. This interpretation was further confirmed by detection a particle with 
energy about 100 EeV 161] The ankle as presented by J. Linsley in 1973 13 
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is shown in Fig. [T] according to data of Volcano- Ranch and Haverah-Park. 
At present beginning of ankle is found at E°^^ ^ 4.5 ± 0.5 EeV according 



to HiRes observations [l7j, at (4.9 ± 0.3) EeV in Telescope Array (TA) [18| 



and at (4.2 ± 0.1) EeV in Auger (PAO) [19|, [20 . 

What makes us think that CR at highest energies are extragalactic and 
that transition occurs at ankle? 

If CRs at highest energies are galactic heavy nuclei up to Iron, anisotropy 
can be not easily observable. However, SNR as the sources cannot provide 
particles with energies as high as 100 EeV, and for all-galactic-CR model 
one needs additional class of sources which are able to accelerate particles to 
much higher energies than SNR. These sources can be hypernovae explosions 
(GRBs) which occur very rarely, e.g. one per million years in our galaxy. 



Such a model was developed in [21|. More recently another galactic model 
[2^ also with GRBs as the sources was studied, with protons and Iron nuclei 
as accelerated particles and assuming diffusion of Iron nuclei in the galactic 
magnetic fields. The calculated energy spectrum explains well the observed 
PAO spectrum. 

To prove that observed particles at highest energies are extragalactic, one 
must know their signatures, and these signature have been found theoretically 



in famous work by K. Greisen [23| and G.T. Zatsepin and V.A. Kuzmin 
24| . who predicted that extragalactic protons propagating through CMB 
radiation obtain the specific steepening of energy spectrum, called Greisen- 
Zatsepin-Kuzmin (GZK) cutoff. 

The highest energy feature, the steepening ('cutoff') of the spectrum is 
found indeed in all three largest detectors, HiRes jisj , Telescope Array (TA) 



261 ] and PAO [27[, though the nature of this cutoff is still questionable. 

The spectrum steepening is quite different for protons and nuclei as pri- 
maries. The energy losses of these particles interacting with CMB and EBL 
(Extragalactic Background Light) photons are shown in Figj2l 

The GZK cutoff is caused by photo-pion production on CMB photons: 

P + 7cmb ^ N + 71 (1) 

The energy Ecq2 ~ 60 EeV (see Fig. [2]) describes this steepening as interaction 
signature. In realistic calculations beginning of the GZK cutoff is lower, 
because at larger z the intersection point Ecq2{z) is shifted by factor (l + z) to 
lower energies. The spectrum must show steepening at lower energy Eqzk ~ 
50 EeV, though this is (weakly) model-dependent quantity. However, the 
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Figure 2: Energy losses of protons (left panel) on CMB and EBL. One can see that 
energy losses on EBL can be neglected. There are two characteristic energies: i?oqi = 
2.4 EeV (intersection of pair-production e^e~ and adiabatic losses) and E^cfi = 61 EeV 
(intersection of pair- and pion- production curves) . It is clear that beginning of the GZK 
cutoff must be about (but below) i?q2. In the left panel the inverse lifetimes of two nuclei, 
Nitrogen and Iron, are shown as function of Lorentz factor F from 28|. The dotted 
curve shows contribution of EBL. The small-dots line shows the case of CMB only, which 
corresponds to steepening at intersection of pair-production curve with adiabatic one . 
Presence of EBL strongly modifies this picture. 



PAO data for beginning of steepening logE'czK = 19-41 ± 0.02 
eV, (or Egzk = 25.7 ± 1.2 EeV) given in Table 1 of 
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with E in 
is below all known 



theoretical predictions, even after taking into account a possible 22% energy 
shift due to systematic energy errors. 

On the other hand HiRes 25[ and TA 26| data agree well with predicted 
beginning of the GZK cutoff and its energy shape. This contradiction reflects 
the different mass composition found in these experiments. In the discussed 
energy region the mass composition of PAO is presented by heavy nuclei and 
the observed spectrum cutoff might be explained as nuclei- destroying cutoff 
(see nuclei life-times in the right panel of FigJ2]). Historically the GZK cutoff 
is related to process p -|- 7cmb — )■ iV + vr. To rename it including there also 
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nuclei-destroying process would be incorrect, because Greisen who shortly 
discussed this process, considered it only for CMB photons as target. In this 
case beginning of the cutoff and its shape is completely different from the 
case of EBL, as one can see it from from right panel of Figj2l 

In contrast to GZK cutoff the nuclei-destroying cutoff has no well pre- 
dicted observational signatures: in particular the energy position and shape 
of the cutoffs are different for different nuclei. 



2. Signatures of proton propagation through CMB 

Propagating through CMB, UHE protons undergo two interactions: photo- 
pion production p-|-7cmb N+n and e+e" production p-|-7cmb p+e^+e~ . 
As a result, the proton spectrum is distorted: due to the first interaction it 
obtains a sharp steepening called GZK cutoff and due to the second one - 
shallow deepening, called dip. Both features depend not only on interactions 
but also on model-dependent quantities, e.g. on modes of propagation (dif- 
fusion or rectilinear propagation), on cosmological evolution of the sources, 
on source separation etc. These model-dependent distortions are especially 
strong for GZK feature. 



The main strategy of this presentation as well as of works t3ul . \3A \32J is 
to distinguish the interaction signatures from the model- dependent ones. 

This is possible to do using the modification factor r](E) js^, some kind 
of theoretical spectrum, in which model-dependent features are suppressed 
or absent. For this aim we use the calculations which involves only one free 
parameter. One naturally expects that the calculated interaction signature 
in terms of modification factor cannot have the agreement with observations 
with good because the observational data include the model-dependent 
features described by many parameters, such as cosmological evolution of the 
sources, source separation etc. One free parameter is not enough to describe 
4-5 different experiments with about 20 energy bins in each. As the next 
step we perform the model-dependent calculations which necessarily include 
many free parameters improving further the agreement. This analysis should 
be done in terms of natural spectrum, J{E) or E^J{E), where the model- 
dependent features are not suppressed. 

However, the Nature has been more kind to us than we expected. The dip, 
in terms of modification factor with one free physical parameter (interaction- 
signature des crip tion) gave very good agreement with observations of four 
experiments [3J]: Yakutsk, AG AS A, HiRes and Telescope Array (see section 
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12. ip . The comparison with PAO data has a different story. Comparison in 
terms of modification factor with PAO observational data from ICRC 2007 



(Merida) had good enough (see 35|). However, with increasing statistics 



the agreement became worse, and comparison of modification factor with 
PAO data 2010 - 2011 has bad x^- It could be that for the difference in 
very small error bars the model-dependent effects are responsible. Indeed, 
as demonstrated in section 12.21 the 22% shift of PAO energy scale and cos- 
mological evolution make the PAO spectrum in terms of E^J{E) compatible 
with the pair-production dip with high accuracy. However, this interpreta- 
tion contradicts the nuclei mass composition measured by PAO in energy 
region of the dip. 

2.1. The dip in terms of modification factor 

In this section the dip will be studied as a signature of UHE proton inter- 
action with CMB, using the modification factor as a tool. The modification 
factor ri{E) is defined as the ratio of proton spectrum Jp{E) calculated with 
all energy losses to the so-called unmodified spectrum J^^^^{E) in which only 
adiabatic energy losses (red-shift) are included: 

r){E) = Jp{E)/J^,,^{E). (2) 

Modification factor is an excellent tool for interaction signature. As one might 
see the interactions enter only numerator and thus they are not suppressed in 
ri{E), while most other phenomena enter both numerator and denominator 
and thus they are suppressed or even cancelled in modification factor, the 
property being especially pronounced for the dip modification factor, which 



according to our calculations [32| depends very weakly on generation index 
7c, and -Emax, on propagation mode, source separation within 1-50 Mpc, 
local source overdensity or deficit etc. The dip modification factor is modified 
strongly by presence of nuclei ( > 15%). Modification factor for GZK feature 
is changing stronger. 

Theoretical modification factors calculated for different source generation 
indices 7^ are presented in Fig. [3l If one includes in the calculation of Jp{E) 
only adiabatic energy losses, then, according to its definition, ri{E) = 1 
(dash-dot line in Fig. [3]). When e"'"e~-production is additionally included, one 
obtains spectrum rj{E) shown in Fig.Oby the curves labeled as ?7ee. With the 
pion photo-production process being also included, the GZK feature (curves 
"total" ) appears. The observable part of the dip extends from the beginning 
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E, eV 

Figure 3: Modification factors for two generation indices 7^ = 2.7 and 2.0. Tire dotted 
curve slrows rj^c when only adiabatic and pair-production energy losses are included. The 
solid and dashed curves include also pion-production losses. 



of the GZK cutoff at £" ^ 40 EeV down to £" ^ 1 EeV, where r] ^ 1. 
It has two flattenings: one at energy E^^ ~ 10 EeV and the other at Eb ~ 
1 EeV. The former automatically produces the ankle (see Fig. H]) and the latter 
provides an intersection of the flat extragalactic spectrum aX E < 1 EeV with 
the steeper Galactic one. We discussed above the theoretical modification 
factor. The observed modification factor, according to definition, is given by 
the ratio of the observed flux Jobs(-E') and unmodified spectrum Junm(-E) oc 
E~^3^ defined up to normalization as: ?7obs oc Johs{,E) / E~^^ ■ Here 7^ is the 
exponent of the generation spectrum Qgen{Eg) oc E~'^^ in terms of initial 
proton energies Eg. Fig. H] shows that both the pair production dip and the 
beginning of the GZK cutoff up to 80 EeV are well confirmed by experimental 
data 3J] of Akeno-AGASA, HiRes, Yakutsk and Telescope Array (TA). 

The comparison of the theoretical dip with observational data includes 
only two free parameters: exponent of the power-law generation spectrum 7^ 
(the best fit corresponds to ■jg = 2.6 — 2.7) and the normalization constant 
to fit the e~^e~-production dip to the measured flux. The number of energy 
bins in the different experiments is 20 — 22. The fit is characterized by 
X^/d.o.f. = 1.0 — 1.2 for AGASA, HiRes and Yakutsk data. This is a very 
good fit for signature (see beginning of this section). For this fit we used 
the modification factor without cosmological evolution of sources. As was 
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Figure 4: The predicted pair-production dip in comparison with Akeno-AGASA, HiRes, 
Yakutsk and Telescope Array data [IJl- AH these experiments confirm the dip behavior 
with good accuracy, including also the data of Fly's Eye [sJl (not presented here). 



explained above, using a model approach with additional three parameters 
describing the cosmological evolution one can further improve the agreement. 
In Fig. m one can see that at E < 0.6 EeV the experimental modification 
factor, as measured by Akeno and HiRes, exceeds the theoretical modification 
factor. Since by definition the modification factor must be less than one, 
this excess signals the appearance of a new component of cosmic rays at 
E < Eij. ^ 0.6 EeV, which can be nothing else but the Galactic cosmic 
rays. This interpretation is confirmed by transition of heavy component 
to the protons in the upper-left panel of Fig. [71 that with good accuracy 
occurs at the same energy. Thus, according to HiRes data the transition 
from extragalactic to Galactic cosmic rays, occurs at energy E^,. ~ 0.6 EeV 
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and is accomplished at i? ~ -E;, ~ 1 EeV (see upper-left panel in Fig. [7] as 
example) . 



2.2. Pair-production dip as energy calibrator 

The energy position of pair-production dip is rigidly fixed by interaction 
with CMB and thus it can serve as energy calibrator for the detectors. 

As we already mentioned, it is difficult to expect that in terms of the 
modification factor the dip described by one free physical parameter can fit 
the observational data with minimum x^- One can shift the observed energy 
bins by the recalibration factor Acai, within the systematic error of observa- 



3l|. We shall refer to this procedure as 'recalibration 



tions, to mmimize 
of energy scale'. We discuss first the dip in PAO spectrum 



29| presented 




Figure 5: Left panel: Comparison of the PAO energy spectrum (filled boxes) with the 
HiRes and TA data fitted by theoretical pair-production dip (solid curve). Right panel: 
Spectra after energy recalibration of the PAO data with Acai = 1-22 (see the text). 



by the filled boxes in the left panel of Fig. [51 As it was already mentioned, 
because of very small statistical error bars it has too large comparison 
with pair-production dip terms of modification factor. We shall use then 
the model-dependent method in terms of E^J{E) including the cosmological 
evolution (1 -|- z)"^ up to ^max as shown in the left panel of Fig. O Using 
two more free parameters m and ^max we can reach better agreement with 
the modified shape of the dip shown the solid curve in Fig. [51 Now we can 
shift the PAO energy bins by factor Acai reaching the minimum x^. For this 
Acai = 1-22 is needed. As a result we obtain picture shown in the right panel 
of Fig. [5l We obtained not only the excellent agreement with the shape of 
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E, eV E, eV 

Figure 6: Left panel: Original fluxes from all detectors (fluxes from HiRes and TA are ap- 
proximately the same). Right panel: Spectra after energy recalibration by pair-production 
dip: Acai ^ 1 for HiRes/TA, Acai = 1.22 for PAO, Acai = 0.75 for AGASA and Acai = 0.625 
for Yakutsk. 



the theoretical pair-production dip (sohd curve) but also the good agree- 
ment with absolute fluxes of HiRes and TA. Note that disagreement with 
GZK cutoff remains for the three energy bins in energy interval 35 — 52 EeV. 
Recalibration with help of pair-production dip for all five detectors (HiRes, 
Telescope Array, PAO, AGASA and Yakutsk) is shown in Fig. [6l Recalibra- 
tion factor Acai = 1 for HiRes/TA is based on the scale factor which correctly 
describes the pair-production dip and GZK cutoff in differential and integral 
{E1/2) spectra. 

2. 3. GZK cutoff in HiRes and Telescope Array data 

The two largest Extensive Air Shower (EAS) detectors, HiRes 25|] and 
Pierre Auger Observatory ji^] have observed a sharp steepening in the UHECR 
spectrum at > (30 — 50) EeV. Both collaborations claimed that the ob- 
served steepening is consistent with the GZK cutoff. But as a matter of fact, 
there is a dramatic conflict between these two results, which still leaves the 
problem open. In this subsection we analyze data of the HiRes which pro- 
vide a strong evidence in favour of the GZK cutoff. These data are supported 



also by the TA data [26l|. The data of PAO will be considered in the next 
subsection. 

To interpret convincingly the spectrum steepening as the GZK cutoff one 
must prove that (z) energy scale of the cutoff and its shape correspond to 
theoretical predictions and (m) the measured mass composition is strongly 
dominated by protons. In HiRes the mass composition is determined from 
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Figure 7: Mass composition and GZK cutoff as measured by tlie HiRes detector. In two 
upper panels (Xmax) (left) and RMS (right) are presented as function of the energy. Both 
agree with a pure proton composition, shown by curves labeled 'proton'. The left- lower 
panel shows differential energy spectrum in terms of the modification factor. One can 
see a good agreement with the predicted shape of the GZK cutoff. The right-lower panel 
shows the quantity Ei /2 in the integral spectrum. This energy, a characteristic of the GZK 
cutoff, is found as £'1/2 = ^O^^'^'^^" ^^ eV in good agreement with theoretical prediction 
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2019-72 gv (see the text). 



(-^max) (E), average depth of atmosphere in g/crn?, where a shower with 
energy E reaches maximum, and RMS(Xmax); which is the width of the dis- 
tribution over Xmax- These values measured by the HiRes are displayed in 
Fig. [71 From the left-upper panel of Fig. [7] one can see that the chemical com- 
position changes from very heavy elements, probably Iron, at ~ 0.1 EeV, 
(data of HiRes-MIA |36|) to protons at E ~ 1 EeV (data of HiRes [ij). 
RMS(Xmax), a very sensitive tool for mass composition, also provides ev- 
idence for a proton-dominated composition at > 1 EeV and up to the 
highest energies (see upper- right panel of Fig. [7]). Differential energy spec- 
trum of the GZK feature in the form of modification factor (left-lower panel) 
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Figure 8: Le/i panel: Auger data 37[ for (Xmax) as function of the energy (left panel) and 
for RMS(Xi„ax), the width of the distribution over Xi„ax, (right panel). The calculated 
values for protons and Iron are given according to QGSJETOl [Sg and QGSJET II [s^ 
models. One can see from the right panel that RMS distribution becomes more narrow 
with increasing energy which implies a progressively heavier mass composition. 



is in a reasonably good agreement with the theoretical prediction, though 
better statistics at higher energies is still needed for a final conclusion. 

The integral energy spectrum of UHE protons, Jp{> -E), has another spe- 



cific characteristic of the GZK cutoff, the energy E1/2 [33[. It is based on 
the observation that the calculated integral spectrum below 50 EeV is well 
approximated by a power-law function: Jp{> E) oc . At high energy 
this spectrum is steepening due to the GZK effect. The energy where this 
steep part of the spectrum equals to the half of its power-law extrapolation, 
Jp(> E) = KE~^, defines the value of -E'i/2- This quantity is found to be 
practically model-independent; it equals to E1/2 = 10^^'^^eV ~ 52.5 EeV 33 . 
Fig. [7] demonstrates how the HiRes collaboration found E1/2 from observa- 



tional data j40|. The ratio of the measured integral spectrum J(> E) and 
the low-energy power-law approximation KE~^ was plotted as a function of 
energy. This ratio is practically constant in the energy interval 0.3 — 40 EeV, 
indicating that the power- law approximation is a good fit, indeed. At higher 
energy the ratio falls down and intersects the horizontal line 0.5 at the energy 
defined as -E'i/2- It results in i?i/2 = 10^^-^^='=°-°'^ eV, in an excellent agreement 
with the predicted value. 

Thus, one may conclude that the HiRes data presented in Fig. [7] indicate 
the proton-dominated chemical composition and the presence of the GZK 
cutoff in both differential and integral spectra. The conclusion about proton 



composition is further supported by the recent TA data [26 
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2.4- PAO data: energy spectrum and mass composition 

In subsection 12.21 we demonstrated that the dip shape, as observed by 
PAO, can agree after recahbration with energy spectra of HiRes/TA and 
other detectors (see Figs. O and E])- The agreement of the PAO and HiRes/TA 
spectra is related to low energy part of the energy spectrum in the right 
panel of Fig. [51 At higher energies statistical uncertainties are too large to 
distinguish between the spectra. 

While the HiRes and TA spectra are compatible with the GZK cutoff, 
the Auger spectrum is not. The steepening in the PAO spectrum starts at 



energy E = 25.7 ± 1.2 EeV |20|, |29l, lower than Egzk ^ 50 EeV, and in 
three successive energy bins in the interval 35 — 52 EeV the PAO flux is 
significantly lower than one predicted for the GZK shape as shown in the 
right panel of Fig. |5l We could not reconcile the PAO cutoff shape with the 
GZK behavior by including in the calculations different generation indices 7^, 
evolution regimes, low acceleration maximum energy -Emax? local overdensity 
of sources etc. 

This disagreement is quite natural for the PAO mass composition data 
which, in contrast to HiRes and TA, show strong dominance of nuclei (see 
Fig. [H]). A steepening in the end of the of the nuclei spectrum, as calculations 
show, is quite different from that of protons (GZK cutoff). The most reliable 
data on mass composition is given by elongation curve Xmax(-E') and espe- 
cially by RMS(Xinax (-£'))• In the Auger data the latter steadily decreases with 
energy and approaches the Iron value at ^ 35 EeV. Low RMS, i.e. small 
fluctuations, is a typical and reliable feature of the heavy nuclei composition. 
These data are further strengthened by other PAO measurements provided by 
surface detectors. They allow to extract two other mass-composition depen- 
dent quantities: the atmospheric depth (X^^x)) where muon-production rate 
reaches maximum, and maximum zenith angle ^9max determined by the sig- 
nal rise-time in surface Cerenkov detectors. Measurements of both quantities 
confirm the heavy mass composition and its dependence on energy obtained 



with the help of (Xmax) (E) and RMS(Xmax), |4l| and [42|, respectively. 

Our further analysis of the Auger spectrum and mass composition is based 
on the following two observations: 

(i) According to the HiRes (Fig. [7]) and PAO (Fig. [8]) data, the observed 
primaries at energy (1 — 3) EeV are predominantly protons or nuclei not 
heavier than Helium. 

These particles cannot be galactic, otherwise, as MC simulations j43| show. 
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galactic anisotropy would be too large. 

(a) The particles at higher energies are extragalactic nuclei with the charge 
number Z increasing with energy. 

This observation is naturally explained by rigidity dependent acceleration in 
the sources Ef^^^ = Zi x E^^^, since at each energy E = ZE^^^ the con- 
tribution of nuclei with smaller Z' < Z vanishes. It was demonstrated in 



4J] that to avoid a proton dominance at the highest energies, one must 



assume that the maximum energy of the accelerated protons is limited, 
£^ax < _ EeV and thus E^^^ ~ 200 EeV. The calculated spectrum 



4J] explains well the PAO energy spectrum, however a problem remains to 



explain simultaneously energy spectrum and mass composition. Somewhat 



similar calculations are performed in j45[ (see Fig. 4 there), and j46 



3. Status of ankle as transition from galactic to extragalactic cos- 
mic rays 

As was discussed in section [T] there are two features in the observed CR 
spectrum, where transition from galactic to extragalactic cosmic rays is ex- 
pected: the second knee and ankle. 

The second knee at energy Egcn ~ (0.4—0.7) EeV as transition feature cor- 
responds well to the Standard Model (SM) of Galactic CRs where maximum 
energy of acceleration corresponds to Iron knee, -Emax ~ (0.08 — 0.1) EeV. 

The ankle is located at E'q ~ (4 — 5) EeV according to all data, and 
as a feature of transition it needs another high-energy galactic component 
beyond SM. On the other hand the ankle appears automatically as the dip 
feature in Akeno-AGASA, HiRes, Yakutsk and TA data (see Fig. H]) and also 
in Auger data (Fig jS] right panel) after recalibration of energies with A = 1.22 
and inclusion of evolution of the sources in theoretical prediction. At this 
comparison the power-law spectra approximation is not needed. Note that 
ankle in this case is a feature produced by interaction of extragalactic protons 
with CMB. We will demonstrate below that interpretation of the ankle as 
transition feature contradicts the observational data of the Auger, HiRes and 
TA detectors. 

The key data used in a proof is given by energy interval (1 - 3) EeV 
where the observed mass composition in HiRes, TA and Auger agrees with 
pure proton composition and cannot be heavier than Helium (for Auger data 
see Fig. [8]). In the ankle model of transition this interval is below the ankle 
at Ea ^ {4: — 5) EeV and thus must have galactic origin. Meanwhile, the 
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Monte-Carlo simulations (see e.g. 47|) and most convincingly the recent one 



43[, show that predicted anisotropy exclude the galactic protons and light 
nuclei at > 1 EeV as observed primary particles. This conclusion is further 
strengthened by recent upper limit on the anisotropy observed by the Auger 
collaboration 



In fact, the conclusion about ruling out the ankle as transition feature 



has been already made in [49( and more recently in [50| from analysis of 
elongation curve Xm»AE)- One can see from right-upper and right-lower 
panels of Fig. 10 in |50(], the strong contradiction of the ankle model with 
the Xinax(-E') data of Auger and HiRes in (1 - 5) EeV energy interval. 

The analysis above proves that transition from galactic to extragalactic 
CRs occurs below 1 EeV, where the only visible spectral feature is the second 
knee. Together with the proton mass composition at (1-3) EeV observed in 
HiRes, TA and Auger experiments it gives one more support to the pair- 
production nature of the observed dip and the dip model developed in [s^, 



31, 32 



An alternative explanation of the ankle observed in PAO at 4 EeV is 
given by transition at this energy from extragalactic protons to extragalactic 
nuclei. However according to PAO observations ji^] the energy of this tran- 
sition using the change in elongation rate and in RMS(E) is given by Et^ = 
2.40lo:72 EeV to be compared with position of ankle Ea = (4.17 ± 0.1) EeV. 
If to assume i^tr ~ Ea the Auger spectrum and mass composition can be 
explained by disappointing model 4^. The transition from galactic CRs 
(Iron) to extragalactic protons occurs at the second knee. Transition from 
extragalactic protons to extragalactic nuclei occurs near the maximum ac- 
celeration energy for protons E^^^ ^ 4 EeV. This energy provides transition 
from protons to heavier nuclei because their maximum energy E'^^^ = ZE^ 
is higher. The observed highest energy steepening can be explained as nuclei 
destroying cutoff. 



max 



4. Conclusions 

The most important result of analysis performed in this paper is the con- 
clusion on interpretation of ankle, observed at ~ 4 EeV in all experiments. 
This feature is usually interpreted as transition from galactic to extragalactic 
CRs. The presented analysis is based on consensus about mass composition 
measured by HiRes, TA and Auger detectors, which show that at energies 



16 



(1-3) EeV (i.e. below the ankle) the primaries are protons or nuclei not heav- 
ier than Helium. Having galactic origin, these particles must show the big 
anisotropy jisj in contrast to low anisotropy measured recently 48 1 by Auger 
detector. In fact this analysis only enhances the similar conclusion reached 



in papers |49| and |50| about ruling out the models with transition at the 
ankle, using the measured elongation curves Xmax(-E')- 

In the light of these results, the transition must occur at lower energies, 
most probably at the second knee from Iron galactic component to proton 
extragalactic component as the dip model predicts 3^ 31, 32|. It strengthens 
further the dip model. 

The dip model is based on the dip cLS cL SI gnature of e~^e pair-production 
in interaction of UHE protons with CMB: p + '~fcmh p + e~ + e~^. The 
characteristic which is sensitive mostly to interaction is given by modifica- 
tion factor r]{E) = Jp{E) / J^nui{E) , where unmodified proton flux Jnnm{E) is 
calculated only with adiabatic energy losses taken into account, while Jp{E) 
is calculated accounting for all energy losses. One may see that interactions 
enter only numerator and thus modification factor is very sensitive to interac- 
tions, while most other phenomena enter both, numerator and denominator, 
and thus they are suppressed or even cancelled in modification factor. The 
theoretical modification factor is presented in Fig. [3] and one can see that it 
practically does not depend on the generation index 7^. Another property of 
modification factor is rj^E) — > 1 at — > 1 EeV. Comparing the theoretical 
and experimental modification factors (see Fig. Hj) one can use only one free 
physical parameter for description about 25 bins in each of four experiments, 
however agreement is quite good. It confirms that observed dip is really 
produced by p -|- 7cmb — )■ p + e~ + e"*". 

One naturally expects that the calculated interaction signature in terms 
of modification factor cannot have the agreement with high-statistics ob- 
servations with very good x^? because the observational data include the 
model-dependent features described by many parameters, such as cosmologi- 
cal evolution of the sources, source separation etc. When statistics increases, 
agreement with modification factor must become worse as it happened with 
Auger data. As the next step we performed the model- dependent calcula- 
tions which necessarily include many free parameters describing all important 
physical effects. This analysis should be done in terms of natural E^J{E) 
spectrum, where the model-dependent features are not suppressed. 

The dip has the fixed energy position and shape, therefore it can serve as 
natural energy calibrator. Shifting the energy scale in each four experiments 
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by factor A, i.e. shifting each energy bin in each experiment a by factor Aq, 
one can reach Xmm fo'^ each experiment. Interestingly enough that as a result 
of this procedure the Auger data show the dip (see right panel of Fig. [5]) and 
fluxes of all experiments coincide, as one can see in Fig. |6]and in Fig. 11 of 



Ref. [3l|. We assume there A = 1 for HiRes experiment, because it has the 
correct scales in the form of beginning of GZK cutoff and -Ei/2- 

The dip is valid for proton-dominated spectrum. Mixing of protons with 
15% of nuclei modifies the dip shape. 

The different mass composition presented as a result of direct measure- 
ments in the HiRes and TA experiments, on one side, and in Auger experi- 
ments, on the other side, is at present the main conflict in UHECR. While 
the data of HiRes/TA show the proton-dominated mass composition at all 
energies above 1 EeV, a few different methods of mass-composition measure- 
ment in Auger experiments demonstrate the nuclei-dominated composition 
starting from energy E 4 EeV and being steadily heavier as energy in- 
creases. Data of all three experiments, HiRes, TA and Auger, coincide in 
narrow energy interval (1-3) EeV, where composition is proton-dominated. 
Each of two possibilities can be viable within two contradicting models. 

The model for description of HiRes and TA data is a proton-dominated 
one with transition from galactic (Iron nuclei) to extragalactic (proton) com- 
ponent at the second knee. The presence of two features, the pair-production 
dip and GZK cutoff, are compulsory for such mass compositions, and these 
features are observed indeed as the shape of the dip and beginning of the 
GZK cutoff, together with value of Ei/2- The dip model describes all these 
features numerically. 

The Auger mass composition implies quite different scenario. The tran- 
sition from galactic to extragalactic CRs occurs at the second knee (-Eskn ~ 
(0.4 — 0.7) EeV) as the only visible spectral feature between the Iron knee 
(galactic feature) and the ankle (extragalactic feature). The ankle at Ea ~ 
4 EeV may be produced by transition from extragalactic proton component 
at < 4 eV to extragalactic heavy nuclei component, though according to 
Auger data 37| this transition occurs at lower energy E 2A EeV. Pro- 
ton component disappears because of low E"™*^^ ~ 4 EeV at acceleration 



44l |. and heavier nuclei progressively dominate at higher energies because 
of the higher acceleration energies E^^^ = ZE^^^, or because of diffusion, 
or both. The highest-e nerg y spectrum suppression, observed by Auger at 
-^'supp = 25.7 ± 1.2 EeV 29| is too low for GZK cutoff. It can be explained 



as nuclei-destroying cutoff. 
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The Auger scenario is characterised by very specific behaviour of X^g^^{E) 
and RMS(E) curves since protons appear in a relatively small energy range 
(1-3) EeV. In this range Xmax and RMS reach the highest values, typical for 
protons, with lower values starting at both ends of these interval, because 
of transition to galactic Iron (lower end) and extragalactic nuclei (upper 
end). Intersection of galactic spectrum (Iron) with extragalactic protons 
corresponds to maximum of RMS. 

All three biggest experiments use the fiuorescent detectors as the main 
tool of measurements. The contradicting results are connected with some 
difference in the data analysis. Close collaboration between Telescope Array 
and Auger teams may resolve this problem in the nearest future. The role 
of new experiments is less significant: already now there two experiments 
which confirm proton-dominated composition and one - nuclei-dominated 
composition. Science is not elections. 
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